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Signal Processing and Systematic Phase Coding for 
Mitigation of Range and Velocity Ambiguities in Doppler 

Weather Radars 

BACKGROUND OF THE INVENTION 

This invention relates, in general, to Doppler weather radar range 
and velocity ambiguity resolution. In particular, (a) to phase coding and 
special processing of uniform PRT sequences for the mitigation of range 
ambiguities as well as, (b) to similar processing to mitigate velocity 
ambiguities and to filtering the ground clutter in staggered PRT sequences. 

Description of the Prior Art 

Range and velocity ambiguities are inherent in a pulse Doppler weather 
radar that operates at centimeter wavelengths. For a pulse repetition time 
(PRT), of T, the unambiguous range is cT/2 where c is the speed of light. 
The unambiguous velocity is given by X/(4J), where X is the wavelength. 
To achieve a large unambiguous range requires a long T, whereas a large 
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unambiguous velocity is inversely proportional to T. Hence the two 
requirements are contradictory. 

The problem is exasperated in surveillance radars at low elevation 
angles because the required pulse repetition time to clear all range 
ambiguities would cause excessive aliasing of Doppler velocities. Current 
operational radars of the NWS (National Weather Service) use two scans at 
low elevations, one at low PRT for a velocity measurement and the other at 
a higher PRT for a reflectivity measurement. Thus, reflectivities are not 
compromised but in the presence of widespread precipitation, velocity 
measurement is usually possible in only one unambiguous range interval 
or annulus of width cT v 12. T v is the PRT for velocity measurement and it 

is usually 4 times shorter than the time T z for reflectivity measurement. 
In extreme cases of wide spread precipitation, echoes from four annuli 
(cT v /2 wide) could be overlaid over each other. 

Methods to mitigate range ambiguities are numerous and one of these 
involves random phase coding. In this method the phases of transmitted 
pulses are coded, purposely, as in the case of klystron transmitters or, 
inadvertently as in the case of magnetron transmitters. The code whitens 
the spectrum of the transmitted (and received) sequences. By recording 
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the coded phase it is possible to reconstruct a coherent pulse sequence for 
echoes from any of the unambiguous range intervals. Thus, reconstruction 
of the stronger echo allows determination of its Doppler spectral moments. 
Subsequent filtering of the strong echo is done and the phases of the 
weaker echo are cohered so that its spectral moments can be computed. 
One of the weaknesses of the method is that filtering of the strong signal 
removes a significant portion of the weak signal. This generates self noise 
in the weaker signal. If only a small part of the strong signal is removed 
its residual will appear as white noise and degrade the spectral moments 
of the weaker signal. 

Another method is the staggered PRT in which the spacing of 
transmitted pulses alternate between Ti and T2 , both of which are chosen 
to be large so that the range ambiguities seldom occur. The unambiguous 
velocity interval is inversely proportional to the difference ITi - T2I. 
However, the variance of the velocity estimates is much larger than what it 
would be if the pulses were transmitted at a uniform PRT equal to 
IT 1 - T2I. Furthermore, no satisfactory procedure to eliminate ground 
clutter from the staggered PRT sequence has been demonstrated. 

In the prior art various methods to resolve range and velocity 
ambiguities in a pulse Doppler radar system have been proposed. For 
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example, U.S. Patent No. 4,730,189 to Siegel et al discloses a Doppler 
radar system with a variable pulse repetition frequency that has a 
coherent integrator to which the reflected pulses are supplied in order to 
avoid deterioration of the indication of a moving target given over-the- 
horizon reception. The coherent integrator subjects the reflected pulses to 
a simplified vectorial addition by the use of filters. 

United States Patent No. 5,150,125 to Hager discloses a Doppler 
radar altimeter designed to resolve the ambiguous range problem 
associated with the use of a pulse repetition interval by transmitting a first 
and second series of pulses where the first series has a pulse repetition 
interval slightly different from the pulse repetition interval of the second 
series. 

United States Patent No. 5,247,303 to Cornelius et al discloses a 
Doppler radar system which reduces ambiguities in measured parameters 
by staggering at least one pulse of each frame of pulses transmitted. The 
staggered pulse creates an isolated pulse doublet pair to provide three or 
more independent and auxiliary Doppler parameter estimates; the 
undisturbed uniform PRT pulses are amenable to ground clutter filtering 
but the doublet pair is not. 

United States Patent No. 5,442,359 to Rubin discloses a method for 
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resolving Doppler frequency shift ambiguities by providing a radiated 
signal that is modulated with a periodic waveform having a plurality of 
pulses within a period with the interpulse intervals in the period being 
unequal. Radar target returns are autocorrelated for a plurality of lags, 
which may include the interpulse intervals and linear combinations 
thereof. 

SUMMARY OF THE INVENTION 

The present invention is a method for resolving range and velocity 
ambiguities in a Doppler weather radar by coding the phases of the 
transmitted pulses in a deterministic manner to generate replicas of the 
spectra from either trips. When the signal from one range interval is made 
coherent, the signal from the other range interval has a multiple split 
spectrum. The spectra have the same shape but are offset from each 
other. The present invention also vastly improves spectral moment 
estimates in staggered PRT sequences and allows efficient removal of the 
ground clutter. 

It is an object of the present invention to provide a new and 
improved method of resolving range and velocity ambiguities in a Doppler 
radar system. 
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It is an object of the present invention to provide a new and 
improved method of resolving range and velocity ambiguities in a Doppler 
radar system by coding the phases of the transmitted pulses of the radar 
in a deterministic manner. 

It is an object of the present invention to provide an improved 
spectral processing procedure to estimate the spectral parameters in 
Doppler weather radar employing staggered PRT transmission to resolve 
range and velocity ambiguity. 

It is an object of the present invention to provide a method to filter 
ground clutter in a Doppler weather radar using staggered PRT 
transmission. 

These and other objects and advantages of the present invention will 
be fully apparent from the following description, when taken in connection 
with the annexed drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 
FIG. 1 is a block diagram of the functioning of a Doppler weather 

radar which uses phase coding. 

FIGS. 2a - 2h is a depiction of a family of codes derived from the 

SZ(8/64) modulation code used in the present invention. 
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FIG. 3a shows the spectrum of the modulation code SZ(8/64). 
FIG. 3b shows a typical Gaussian weather like spectrum obtained 
from simulations using a 512 points discrete Fourier transform. The 
parameters indicated in the figure are, unambiguous velocity v a = 32 ms- 1 , 
mean Doppler velocity v= -12 m s- 1 and the spectrum width is w = I ms- 1 

FIG. 3c shows a spectrum of the signal corresponding to 
FIG. 3b but after modulation with the SZ(8/64) code in FIG. 3a. 

FIG. 3d shows the same view as in FIG. 3a but the weather like 
simulated spectrum has a width of 6 ms- 1 - 

FIG. 3e shows the same view as in FIG. 3c but the modulated signal 
has a spectrum as in FIG. 3d. 

FIG. 4 shows an azimuthal sector of a weather scanning radar. The 
unambiguous range is r a and overlaid echoes originate at r and r + r a . 

FIG. 5 shows a flow chart of significant steps in the procedure for 
recovery of the spectral moments of overlaid echoes. 

FIGS. 6a - 6f show an illustration in the spectrum domain of steps 
that lead to the recovery of the weaker echo's mean velocity and spectrum 
width. 

FIG. 7 shows a proposed sample overlap scheme for obtaining a 64 
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length sequence. 

FIG. 8a shows the standard error in velocity estimates for the 

procedure that uses the SZ code. 

FIG. 8b shows the standard error in velocity estimates for the 
procedure that uses the random phase code. 

Figs. 9a-f show the spectral processing steps in the staggered PRT 
transmission scheme for ambiguity resolution. 
i Fig. IO shows the improvement in the standard error in the velocity 
£ estimate using the present spectral processing method over the method 
4 using the ratio of autocorrelations. 

D Fig. 11 shows the effectiveness of the present clutter filtering 

i procedure in estimating the velocity and resolving the velocity ambiguity. 

DESCRIPTION OF THE PREFERRED EMBODIMENT 

In the proposed implementation the phases of the transmitted radar 
pulses are coded in a deterministic manner to generate replicas of the 
spectra from either trips. That is, when the signal from one range interval 
is made coherent, the signal from other range intervals have multiple split 
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spectra. The spectra have the same shape but are offset from each other. 
Filtering of the strong signal (in the Fourier transform domain) is much 
less detrimental to the desired weak signal than in the case of random 
phase coding. This is because the essential information about the signal is 
contained in the spectral replicas, therefore, as long as two of the replicas 
are preserved during the filtering process, the self noise would be small, 
in this manner it is possible to eliminate a larger portion of the strong 
overlaying signal and achieve a much smaller degradation of the weak 
signal than is possible with the random phase method. That is, the weaker 
signal can tolerate a much wider notch width of the filter. 

An embodiment of the invention is depicted in Fig. 1. This figure is 
applicable to any phase coding scheme, including random phase coding. 
Our contribution is the addition of a digitally controlled phase shifter, the 
code, and the method to process the overlaid signals, 
a) Code 

The code is derived from Chu's code (1972) . Chu suggested a code for 
which the code length M is not divisible by n so that its autocorrelation 
function has a peak at zero lag (and multiples of M) and is zero otherwise. 
Chu's code is cyclic and can be expressed as: 
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(j) k = I17lk 2 / M; where k=0, 1, 2, ... 

In our case M is divisible by n and there are M/n autocorrelation 
peaks within the code length M. This code (|) k is referred to as the 

modulation code. For example if the phase shifter is encoded with a 
suitable code (to be defined shortly) and the first trip echo is phase 
decoded accordingly (i.e. cohered), the second trip echo will exhibit a phase 

sequence (|>k given by the modulation code. After phase decoding the 
phases §\ >m of the ith trip return are given by: 







k ; for the first trip signal (coherent, no modulation); 


<|>2,k 




" *P k ; for the second trip signal; 


<t>3,k 




" ^ k ; for the third trip signal; 


<t>4,k 


= k-3 


' *P k ; for the fourth trip signal. 



The switching code 4* k of the phase shifter is obtained from the sums 
of the sequence (j)k. 
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k 

¥ k = . 2, nilm2/M; k=0, 1, 2, ... M-l 
m=0 

Thus, the proposed switching code is a phase sequence which is a 

summation of the Chu type code. Furthermore, we have generalized the 
class of codes that can be used as a switching sequence which is the 
following code: 

k 2q 

= . X (n7C/M) 2 (m+p) 2 ; q=0, 1, 2, ... (M /2n-l), 
m=0 p=o k=0. 1, 2,...L-1. 

We refer to these as the SZ(n/M) codes. Also for more than two 
overlaid echoes a subset of codes would work as well (e.g., code n=8,M=64). 
For the purpose of ambiguity resolution the processing sequence length L 
should be an integer multiple of the code length M. 

To generate the switching code one starts with a modulation code 

<|) kte) as follows, where q is an integer. 
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2 q 

(|) k (q) = (I171/M) X (k+p)2 

p = 0 

Furthermore, each permissible q yields a different code. Addition of 

a constant to(|)k (q) shifts the modulation code and thus generates a 

switching code with different phases. This, however, does not change the 

code property or its performance. A negative of <|M q) is also a code. The 

modulation code can be conveniently viewed by plotting its phases on a 
unit circle. For the SZ(8/64) this is illustrated in Figs. 2a - 2h. From any 
one code, other codes can be generated by taking a complex conjugate or 
by taking a mirror image of the code with respect to any line bisecting the 
circle. Because these two operations are independent, together they yield 
4 codes. All these codes have similar characteristics as far as the present 
scheme of decoding is concerned. 

The spectrum of the modulation code of the SZ(8/64) code is plotted 
in Fig. 3a. In Fig. 3b the spectrum of a signal is shown that when 
modulated with the code produces the spectrum as shown in Fig. 3c. 
Eight spectral replicas are clearly displayed in the case of relatively 
narrow spectra. Modulation of signals that have wider spectra causes 
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significant spectral overlap as can be seen by examining Fig. 3d and Fig. 
3e. 

b) Retrieval of velocities 

Processing to retrieve spectral moments builds on judicious 
combinations of classical techniques and on the introduction of magnitude 
domain spectral deconvolution. In what follows two unambiguous range 
intervals are considered and it is assumed that the second trip echo is 
weaker. A sector scan shown in Fig. 4 illustrates the two range intervals 
and two locations from which echoes arrive at the receiver at the same 
time. The notation for spectral moments of the two echoes are as follows. 
The power, mean Doppler velocity and spectrum width of the first and 
second trip echoes are pi, vi, wi, and p 2 , v 2 , w 2 , respectively. The 
algorithm to retrieve the spectral moments of the first and second trip 
echoes is described below. It is meant for use by a fast dedicated signal 
processor and it assumes that the transmitted pulses have been encoded 

with an SZ type code. 

The procedure is applied to the digital samples (complex in phase 
and quadrature phase if done at base band) from each range location (i.e., 
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range gate). It is as follows: 

1. For each range gate take the time series Eik ; k=l,2, .... M. 

I The phase switching sequence is H'r; SZ(n/M) code. 

2. Cohere the 1st trip signal. 

•E|(k )= E|(k ) exp {j^k} 

•1st trip is coherent; 2nd trip is phase coded by a sequence 
(|) k = rmK2/M ; k=0, 1, 2 M- 1. 

3. Multiply by von Hann window weights or other taper, hk 

•E| (k) = E| (k) h k 

4. Filter the ground clutter. 

•Ei (k) = filtered (Ei (k)). 

5. Cohere the second trip . 

•E 2 W= E| k exp {j(|)k} 

6. Autocovariance process Ei and E2 to get pi, vi, wi, wf and P2, v 2 , w 2 , w 2 ' 



Page 14 



Docket No. DOC-9801 

(for the computation of wi \ w 2 ' use Eq. 6.27 of Doviak and Zrnic, 1993, 
and for the computation of wi , w 2 use Eq. 6.32 of Doviak and Zrnic, 1993) 
For computation of wi ' and w 2 ' use equation 6.27 of Doviak and 

Zrnic (1993) which is: 
X 

w = I 1nU_J 1 1/2 

2nT s ^2 R(1) 

where w is wi* or W2\ p is the power estimate of Ei or E2, and R(l) is the 
autocovarience estimate at lag 1 of the signal Ei or E 2 , and X is the radai 
wavelength and T s is the pulse repetition time. 

For the computation of wi and W2 use equation 6.32 of Doviak and 
Zrnic (1993) which is: 

X 

w = I 1n I Rm I I 1/2 

27tT s VT R(2) 

where R(2) is the autocovariance estimate at lag 2. 

7. Compute wi'/w2' ratio. 

•if wi7w2'> 1: trip = 2 (second trip is stronger) - process E 2 . 
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•if wi '/W2 1 < 1: trip= 1 (first trip is stronger) - process E\ . 

8. If trip=2, interchange Ei & E 2 , and all the parameters in step number 6. 

•With this interchange, Ei is the time series with stronger signal 
coherent. 

• Recover p2, V2 and w 2 of the weaker signal. 

[Note: The processing steps 9 to 17 are same for the two cases in step 7 
with Ei replaced by E2. This is accomplished by step 8 and the trip 
numbers are restored in step 18.] 

9. Compute spectrum of Ei. 

•Sr = DFT[Ei] 

10. Set to zero notch n w M coefficients centered on vj; the notched 

spectrum is denoted by Fig. 6d. 

Note: (a) n w is not to exceed the maximum permissible value, (l-2n/M). 
(b) for SZ(8/64) & SZ( 12/64) optimum notch filter center 
location to be computed if trip = I (i.e. 1st trip stronger) and 
ground clutter filter is applied. 
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11. Compute mean power p from the remaining coefficients. 
Multiply p by I/(l-n w ) to get mean power P2. 

12. Compute power ratio pr = 10 log io(Pi/P2> dB. 

13. If pr< 25 dB, correct error in pi estimate. 

# Pl '= Pi " P2 

•compute corrected power ratio pr = pi'/p2 

14. Cohere the weaker signal in Si as follows, 
•d = IDFT [Si] 

•if trip = 1, e 2 = ei exp {-j(|)k} 
•if trip = 2, e 2 = ei exp {j(|)k} 

15. Compute autocorrelation R(l) for e 2 , and compute mean velocity, v 2 

16. Magnitude de-convolution, (for SZ(8/64) and SZ( 16/64) only) 
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•compute magnitude spectrum, S2 - I DFT 1 ^) I 
♦Multiply by the de-convolution matrix, S2 = D S2 1 
[The de-convolution matrix, D, is pre-computed and supplied to the 
processor] D is obtained by replacing the elements in the 
convolution matrix with their magnitudes and then inverting the matrix. 

17. Compute autocorrelation R(l) for S2, and compute width, W2. 

18. If trip = 2, inter-change parameters (pi,vi,wi) and (P2 ? V2,W2> 

19. Store or display the 1st and 2nd trip parameters then process signals 
from the next range gate. 

A flow chart of the processing procedure is shown in Fig. 5. 

An example from simulated time series data is shown in Fig. 6 which 
illustrates a few essential steps in the recovery of the velocity of the 
weaker signal. 

The spectrum of the first trip echo in the absence of the other 
overlaid echoes is in Fig. 6a. Similarly the spectrum of a much weaker 
second trip echo with no overlaid signals is in Fig. 6b; the scales on the 
abscissa differ by a factor of 100. If these two signals are overlaid the 
spectrum looks as in Fig. 6c where the weaker signal is not discernible. 
Centering the notch filter on the first trip spectrum and thus eliminating 
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75% of the spectral coefficients results in a spectrum shown in Fig. 6d. 
Cohering the second trip signal produces a wider spectrum than the 
original (compare Fig. 6b with Fig. 6e). The cohering process consists of 
inverse DFT and multiplication by the shifted complex conjugate of the 
modulation code. The spectrum of the cohered signal is shown in Fig. 6e. 
Deconvolution in the magnitude domain of the spectrum in Fig. 6e results 
in a faithful reproduction of the original spectrum of the weaker (second 
trip) signal (compare Fig. 6d with Fig. 6b). 

The procedure owes its success to the code, as expressed above. 
Furthermore there are several crucial original steps that make the 
difference between success and failure. 

First, there is the notch filter which has a maximum normalized 
width of (l-2n/M). Smaller widths remove less of the undesirable signal 
and larger widths preclude reconstruction of the weaker signal. 

Second there is the deconvolution procedure used for obtaining the 
spectrum width. It is performed in the magnitude domain. In this 
procedure the deconvolution matrix is synthesized from the magnitudes of 
the spectral coefficients of the code and is stored in the memory of the 
processor. This is crucial for obtaining the spectrum width of the weaker 
signal. 
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The third element deals with the code and processing of overlapped 
signals. Cyclic codes with lengths that are powers of two are very 
convenient because such sequences are amenable to efficient Fourier 
analysis. Several of our codes have the shift invariant property. This is 
desirable if one needs to preserve azimuthal resolution. Thus weighting by 
a window function and overlapping the signals for processing is natural for 

this code (see Fig. 7). 

Fig. 7 shows the positions of 128 contiguous time samples. The first 
64 are weighted by the Von Hann weights and then processed according to 
the outlined algorithm. Some of our codes allow the processing for the 
next estimate to start before the sequence for processing the current 
estimate has ended. For example, in Fig. 7 new estimates are generated 
every 52 samples but the code length is 64 samples; new samples can be 
generated at other than 52 spacings as well. This is significant because the 
azimuthal resolution which depends on the number of processed samples 
needs to be preserved. 

Figs. 8a and 8 b give a comparison of performance of the present code 
and the previously used random phase code. The procedures for 
retrieving the velocity using these codes are the same except for the notch 
filter width which is larger for the SZ code and the deconvolution for 
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estimating hej spectrum width which is not applicable to the random phase 
code. These figures were obtained from simulations of overlaid signals. 
The signals were encoded with the SZ(8/64) code and subjected to 
retrieval according to the algorithm The signal sequence is composed of 64 
samples. The width of the notch filter is 75% for the SZ code and 50% for 
the random phase code. The spectrum width of the second trip signal is 
4m/s. On the ordinate is the Doppler spectrum width wi of the stronger 
(first trip) signal. On the abscissa is the power ratio of the stronger to 
weaker signal. The 2.5 dB by 0.5 m/s rectangles are shaded to represent 
the average standard deviation of velocity estimates V2 of the weaker 
(second trip signal). The results are valid for a typical wavelength of 10cm 
(pulse repetition time of 0.78 micro seconds) or for a 5 cm wavelength but 
with a two times shorted pulse repetition time. At other wavelengths the 
graphs are valid for pulse repetition times equal to 0.78 micro seconds 
times one tenth of a wavelength (in cm). 

The magnitude domain deconvolution procedure, explained in the 
context of uniform PRT transmission and SZ phase coded Doppler weather 
radar signals, in conjunction with the spectral domain clutter filtering , is 
applied to the signals from a Doppler weather radar using staggered PRT 
transmission to resolve the range and velocity ambiguity. 
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In the proposed implementation, the transmitted pulses are staggered 
in time, i.e., the PRT is changed to Ti and T2 alternately. This scheme 
provides an unambiguous velocity 

1 

V a >)=±'?l/{4(T2-Ti)} 

/ 

and an unambiguous range 
r a = cTi/2. 

The major problem which has prevented implementation of the 
staggered PRT scheme in practical radars is the clutter filtering. All the 
available methods of clutter filtering produce spurious rejection bands in 
the filter response which prevents estimation of the velocity in these 
frequency bands. 

The proposed spectral domain clutter filtering procedure eliminates this 
problem of loss of velocity information. The procedure also uses the 
magnitude domain deconvolution to reconstruct the signal spectrum, which 
allows estimation of the spectral parameters with much 
lower standard error than was possible before. 
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Figs. 9a-9f illustrate the processing steps for signal samples obtained 
in the staggered PRT mode. Fig, 9a shows a simulated signal spectrum 
which is sought to be recovered. In Fig. 9b is the clutter spectrum. The 
clutter power is 40 dB above the signal power. The uniform PRT sequence 
with the spectrum in Fig. 9c is obtained from the staggered PRT sequence 
by inserting zeros for the missing samples. This spectrum consists of 
weather signal and clutter, but the signal is not visible because it is 40 dB 
below the clutter. Note that the clutter power spreads over five bins. Fig. 
9d is the recovered signal after clutter is filtered. The original spectrum is 
reconstructed using the magnitude domain deconvolution as shown in Fig. 
9e. There are some residuals that remain because of the overlap of the 
signal and clutter spectral coefficients. These are deleted after the mean 
velocity is estimated (Fig. 9f). 

The clutter filtering and the spectrum reconstruction procedure is as 
follows: 

1. The received stagger PRT sequence is converted to a uniform 
sequence by inserting zeros for missing samples. Then suitable window 
weights are applied. The discrete Fourier transform of this sequence is 
expressed as a convolution of the original signal (plus clutter) spectrum 
and the spectrum of the code. The code is a sequence of ones and zeros; 
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ones represent the staggered samples and zeros represent the missing 
samples. For example, the code is 1010010100 ... etc. for T1/T2 =2/3. In 
matrix notation this is written as 

[V] = [C v ] [E] 

where [E] is the column vector of the spectral coefficients of the original 
signal (plus clutter), [C v ] which is the convolution matrix obtained from the 
code spectrum, and [V] is the column vector of the spectral coefficients of 
the derived uniform PRT sequence (i.e., inserted with zeros and 
windowed). 

If N is the number of samples in the zero padded sequence, the 
convolution matrix, [C v ], will have only p non-zero coefficients in each 
row(or column), determined by the code periodicity 
(e.g., p=5 for T1/T2 =2/3). 

2. In the proposed method spectral coefficients Vj ; i=l,2.... N, are cut 
into p row vectors, each having M=N/p coefficients, and are written as a 
(pxM) matrix, [V r ]. 

3. The non-zero elements of the code spectrum are arranged as a 
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column matrix, [Ci], of size (px 1), from which the (pxp) clutter filter matrix, 
[ChL is obtained as 

[Cfil = [Ci] [Ci]t. 

The superscript " t" represents complex conjugate transpose of the matrix. 
Similarly, column matrix [C 2 ] is obtained by shifting the elements of [Ci] by 
one; i.e., take the first element and place it at the last position. The second 
clutter filter matrix is obtained as 
[C f2 ] = [C 2 ] [C 2 ] 1 • 

4. Select a width for the clutter filter in terms of the number of 
coefficients in the spectrum. The clutter is confined to the first few and 
the last few columns of the matrix [V r ]. Form two diagonal matrices, [In] 
and [If2] whose diagonal consists of ones corresponding to the columns 
where the clutter is to be filtered and zeros for the rest. The diagonal of 
the first matrix has first few ones, and the diagonal of the second matrix 
has the last few ones and the rest zeros. 

5. Compute the clutter filtered signal spectrum matrix, [V f] , using the 
equation 
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[Vfl = [V r ] - [Cfi] [V r ] [I f i] - [C f2 ] [V r ] [I f2 ] 



This equation subtracts the clutter from the spectrum and retains only the 
weather signal part. This clutter filter does not bias the velocity and has no 
spurious rejection bands. 

6. Rearrange the elements of [Vf] into a column matrix and then carry 
out the magnitude domain deconvolution to reconstruct the original signal 
spectrum. Mathematically, this is written as 



[IEI] = [IC V I]- 1 [IVfl] 

Here, the deconvolution matrix is the inverse of the magnitude convolution 
matrix. 

7. Compute the mean velocity from the magnitude spectrum [ IEI ] either 
using the phase of the autocorrelation or the first moment of the spectrum. 

8. The residual spectral coefficients outside (2N/p -q) coefficients (q is 
the clutter filter width) centered on the mean velocity are deleted and the 
mean power and spectrum width are estimated using standard formulae (q 
is the width of the clutter filter in terms of the spectral coefficients). The 
residuals do not produce appreciable bias error in the velocity. 
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The proposed signal processing method on staggered PRT sequences 
achieves a much lower standard error in the estimates than available from 
existing techniques. Fig. 10 gives comparison of the standard errors in the 
velocity estimates obtained with the present method and the previously 
reported method. Curve (a) is for the ideal situation where a complete 
sequence is available, and is given as a reference. 

The velocity recovery by the present method in the presence of ground 
clutter (40 dB clutter- to-signal ratio) is depicted in Fig. 11. The velocity 
estimates are nearly bias free and there are no spurious rejection bands. 
The large errors at the extreme ends are due to aliasing. 

Although the method for resolving range and velocity ambiguities in 
a Doppler weather radar and the method of using the same according to 
the present invention has been described in the foregoing specification 
with considerable details, it is to be understood that modifications may be 
made to the invention which do not exceed the scope of the appended 
claims and modified forms of the present invention done by others skilled 
in the art to which the invention pertains will be considered 
infringements of this invention when those modified forms fall within the 
claimed scope of this invention. 
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Claim 1. A method for resolving range ambiguities and separating 
overlaid signals in a Doppler weather radar, wherein said method 
comprises: 

phase coding transmitted pulses from said radar with a special code, 

and 

phase decoding of said coded pulses to make echoes form a specified 
range ambiguous interval coherent, and 

said decoding splits the spectrum from overlaid echoes from said 
range interval, 

said spectra have the same shape but are offset from each other. 

Claim 2. The method for resolving range ambiguities as claimed in 
claim 1, wherein a phase *F k of the k th transmitted pulse is determined by 



k 2q 

¥1^= - X (nn/M) I (m + p) 2 ;q=0, 1, 2, ... (M /2n-l), 
m=0 p = o k=0. 1, 2,...L-1. 
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Claim 3. The method for resolving range ambiguities and separating 
overlaid signals in a Doppler weather radar, wherein said method 
comprises: 

transmitting a phase coded sequence, 

for each range gate take the time series sequence, 

cohere a first trip echo, 

multiply by a window weight, 

filter the ground clutter, 

cohere a second trip echo, 

autocovariance process said first and second trip echoes, 

compute the ratio of estimated widths of the widths of said first and 
second trip echoes, 

if said ratio is less than unity, said first trip echo is stronger than 
said second trip echo, 

if said first trip echo is stronger, process the time series samples in 
which the first trip echo is coherent, 

if said second trip echo is stronger than said first trip echo, process 
the time series samples in which the second trip echo is coherent, 

recover the power, mean Doppler velocity and the spectrum width of 
the weaker trip echo. 
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compute the spectrum of the stronger trip echo, 

notch coefficients centered on the mean Doppler velocity of the 

stronger trip echo, 

compute the mean power of the weaker trip echo, 

compute the mean power ratio of the stronger to weaker trip echo, 

if the power ratio is less than 25dB, correct the error in the mean 

power of the stronger trip echo, 

compute the corrected mean power of the stronger trip echo, 
cohere the weaker trip in the spectrum, 

compute autocorrelation for the cohered weaker trip echo, and 
compute the mean velocity of the weaker trip echo, 
compute the magnitude spectrum for the weaker trip echo, 
multiply the magnitude spectrum for the weaker trip echo by the 

de-convolution matrix, 

compute autocorrelation for the weaker trip echo, 

compute the spectrum width for the weaker trip echo, 

display the power, mean Doppler velocity and spectrum width of the 

stronger trip echo and the weaker trip echo, 

output the power, mean Doppler velocity and spectrum width of the 

stronger and weaker trip echoes and proceed to the next range gate. 
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Claim 4. The method as claimed in claim 3, wherein the spectrum 
width of the weaker trip echo is an estimation, and 

said estimation consists of magnitude domain spectral deconvolution, 

and 

said deconvoluted spectra are processed to retrieve the spectrum 
width. 

Claim 5. A method to filter the ground clutter from the signal 
samples of a Doppler weather radar using a staggered PRT transmission 
scheme for velocity ambiguity resolution, wherein said method consists of: 

transmitting pulses with pulse repetition times Tl and T2 
alternately, 

inserting said signal samples with zeros to make it a uniformly 
sampled sequence, and 

then multiplying by a suitable window function, 

rearranging the spectral coefficients of a uniform sample sequence 
into a matrix [V r ] , 

computing a clutter filter matrices [Cfi] and [Cf2], based on a code 
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sequence for the T1/T2 ratio, 

determining the clutter filter width in terms of the number of 
spectral coefficients, 

multiply the column vectors of said matrix [V r ] containing the clutter 
power by said clutter matrix and subtracting the result from said column 
vectors. 

Claim 6. A magnitude domain deconvolution procedure to 
reconstruct the signal spectrum from signal samples of a Doppler weather 
radar using a staggered PRT transmission scheme for velocity ambiguity 
resolution, wherein said method comprises: 

rearranging a matrix [V r ] into a column matrix and then carrying out 
a magnitude domain deconvolution, and wherein 

said magnitude domain deconvolution consists of multiplying a signal 
spectrum magnitude [ IVI ] by the inverse of the magnitude convolution 
matrix, 

computing the mean velocity from the deconvolved spectrum, 
deleting the residual coefficients and then estimating the spectrum 
width and mean power. 
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ABSTRACT 

A method for resolving range ambiguities and separating overlaid 
signals in a Doppler weather radar is proposed. For uniform PRT 
transmission, it consists of a special deterministic code for phases of the 
transmitted pulsed, and associated decoding and processing of return 
signals. In the decoding process when the signal from one range interval is 
made coherent, the signal from the other range interval has a multiple 
split spectrum. The multiple spectra have the same shape but are offset 
from each other. Processing steps to separate the overlaid signals and a 
procedure to estimate the spectral moments are given. One crucial aspect 
in this procedure is the magnitude domain deconvolution. The magnitude 
domain deconvolution is also applied to a staggered PRT transmission 
scheme whereby it enables the estimation of spectral parameters with 
much lower standard error than the known methods. Moreover, the 
magnitude domain deconvolution combined with special spectral 
processing solves the outstanding problem of ground clutter filtering in the 
staggered PRT sample sequence. 
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